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ABSTRACT 

Unlike the random radial orientation distribution of field elliptical galaxies, galaxies in a clus¬ 
ter are expected to point preferentially towards the center of the cluster, as a result of the 
cluster’s tidal force on its member galaxies. In this work an analytic model is formulated to 
simulate this effect. The deformation time scale of a galaxy in a cluster is usually much shorter 
than the time scale of change of the tidal force; the dynamical process of the tidal interaction 
within the galaxy can thus be ignored. An equilibrium shape of a galaxy is then assumed 
to be the surface of equipotential, which is the sum of the self-gravitational potential of the 
galaxy and the tidal potential of the cluster at this location. We use a Monte-Carlo method 
to calculate the radial orientation distribution of these galaxies, by assuming the NFW mass 
prohle of the cluster and the initial ellipticity of held galaxies. The radial angles show a sin¬ 
gle peak distribution centered at zero. The Monte-Carlo simulations also show that a shift of 
the reference center from the real cluster center weakens the anisotropy of the radial angle 
distribution. Therefore, the expected radial alignment cannot be revealed if the distribution of 
spatial position angle is used instead of that of radial angle. The observed radial orientations 
of elliptical galaxies in cluster Abell 2744 are consistent with the simulated distribution. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 2744 - galaxies: 
kinematics and dynamics - galaxies: structure. 


1 INTRODUCTION 

The total mass of a cluster of galaxies is dominated by dark matter, 
which couples only through gravity with ordinary matter. There¬ 
fore, to map the density profile of a cluster, one needs to find 
how luminous matter in a cluster is related to its gravitational po¬ 
tential. Each way of doing so leads to a method of cluster mass 
mapping. Assuming some sort of dynamical balance established, 
one can model how the spatial distribution and velocity disper¬ 
sion of galaxies respond to its gravitational potenti al and then de¬ 
velop a technique to determ ine the mass distribution dCarlberg et al.l 
Il997l : lfent & Gumilll982h . Hydrostatic equilibrium supposed, the 
density and temperature profiles of hot gas are related to the 
shape of the gr avitational potential well jpabricant et al.l Il98d : 
IXue & \^l2002h . Finding the manner of light ra ys bent by gravity 
gives birth to the method of gravitat ional lensing l lKaiser & Squired 
I 1993 I : lUmetsu & Broadhu^ l2008h . Although many approaches 
have been made, the results from those are not always in good 
agreement dWu et al.ll998l : lBMtelmann & SteinmetJl996h . There¬ 
fore, it is still useful to find other independent probes. 

An extended object in a non-uniform gravitational field feels a 


* E-mail: rongyu@ihep.ac.cn 


distorting tidal force. The tidal force from the mass of a cluster can 
change the shape of its galaxies in a way that the radial orientation 
angles (see the upper panel of Fig. [T] for an illustration) of ellipti¬ 
cal galaxies in a cluster show an anisotropic distribution toward the 
cluster center different from the isotropic distribution of field galax¬ 
ies; this phenomenon is referred to as radial alignment (RA). If RA 
really exists and is caused by the tidal effect of the cluster mass, 
then a new link between an observable, i.e. the radial angle distri¬ 
bution (RAD) of elliptical galaxies, and the gravitational potential 
of the cluster can be found. 


The idea that the cl uster tidal force can cause RA was first 
proposed by Thompson jThompsonll 19761) . in order to explain “a 
possible indication that the galaxies are preferentially aligned along 
the radius vector to the c enter of the cluster” for the Coma cluster 
jHawlev & Peeble^ll975h . In recent years, more observational evi¬ 
dence has been found to support this idea fP ereira & Kuhnll^OSl : 
lAgustsso n & Brai nerd 20(361 : iFaltenbacher et alj 2007h . Theoreti- 
callv. lCiotti & Duttal(ll994l) studied the dynamical process of shape 
changing of an elliptical galaxy under the cluster’s tidal field using 
N-body simulations, and showed an alignment tend ency between 
the ga laxies’ major axes with the radial direction. IPereira et al.l 
( l2008h found strong RA of dark matter substruc t ures i n a cluster 
halo in N-body simulations. lUsami & Fuiimot3 ( Il997h studied a 
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gaseous ellipsoid of uniform density orbiting in the logarithmic po¬ 
tential, and discussed the application of their result to the alignment 
of galaxies in rich clusters. They concluded that inside a critical 
radius, galaxies would be disrupted by the tidal force of the clus¬ 
ter, whereas beyond the critical radius the major axes of elliptical 
galaxies are locked in radial direction. None of the previous works 
take the intrinsic ellipticity distribu tion of galaxies into consider¬ 
ation. In lUsami & Fuiimot3 l ll997ti all elliptical galaxies’ intrin¬ 
sic shape were set spherical, unlikely in real cases. In our work, 
the ellipticity distribution of field elliptical galaxies, which can be 
thought unaffected by environment, is used as the intrinsic elliptic¬ 
ity distribution. 

This work aims to formulating an analytic model to simulate 
RA quantitatively. Section 1 is to verify qualitatively that the shape 
changing timescale for most galaxies in a cluster is short enough 
compared with the Keplerian motion period so that we can ignore 
the dynamical process of tidal distortion and assume an equilibrium 
shape of the galaxy. In section 2 we calculate the equilibrium shape 
of an elliptical galaxy analytically. Then the projection effects are 
included and RAD is generated by a Monte-Carlo method. Finally 
the observed RADs of the clusters Abell 2744 are compared. 


the galaxy. 


^grav “ fe ^dm 

^ GMq ^ GMdn, 

" ^ ( 6 ) 

GnpoRo + 


where fe and are the gravity acting on the unit mass from the 
luminous and dark matter at the effective radius of a galaxy; Mq 
and Rq are the mass and effective radius of the luminous matter 
in the galaxy. The luminous matter is assumed homogeneous, with 
density po ~ 10 ®Mq. Tfdm is the enclosed mass of the dark matter 
within the effective radius. In this work, we only consider the tidal 
distorsion within the effective ra dius. The density profil e of the dark 
matter in an elliptical galaxy is ( iRobles & Matosir2012h 
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2 QUALITATIVE DISCUSSION 


Consider a galaxy in a cluster, and the vector from the galaxy center 
to the cluster center is L. The enclosed mass within radius L is 
denoted by M(L) and its potential field is The tidal force acting 
on a mass element of the galaxy is approximately 


d 3<hc 
dfj dr, 


(1) 


where r and f are the vectors from the mass element to the galaxy 
center and to the center of the cluster, respectively, i = 1,2,3 denote 
the three cartesian coordinate components of a vector. The Einstein 
summation convention is used in this work. We define a tidal po¬ 
tential corresponding to the tidal force, which is 


— 



( 2 ) 


The mass distribution of the cluster is treated as spherically 
symmetric, so its gravity acting on a mass element can be written 
as 


a. GM(f) ^ 

®c --- ^ ®m(^)5 
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where M(f) denotes the enclosed mass in radius r, and Omfr) de¬ 
notes the gravit y potential caused by the m ass gradient. Here since 
-V®c = teinney & Tremainal987h . we obtain 


G • VM(f) 

\f\ 


- VOMCr) = 0. 


(4) 


Substituting Equations 0 and a into Equation m , and neglecting 
the variation of M(r) in small scale, we obtain 


(fix = 


GM(L) 1 


(5) 


where I, is the component of r along L, and Z, is the perpendicular 
component, which are illustrated in Fig.[T] Therefore the tidal force 
acting on an unit mass is obtained by fx = -V®x. 

An unit mass element is also subjected to the self-gravity from 


where p^ ~ O.lA/Q/kpc® is the central density, and rc ~ 5 kpc is the 
scale radius. Therefore the dark mass enclosed within r is = 
4np,.rl(j- - arctan(^)). At the effective radius Rq ~ 10 kpc, == 
GnpaRa ~ 10®, fdm - GM^^IRq ~ 10®, where ttG = 1 is used. 
Therefore fdm « A, fgrav - fe, implying that the effect of the dark 
matter in an elliptical galaxy can b e neglected within the effective 
radius . This has been comfirmed bv ide Paolis. Ingrosso & Strafellal 
( Il995h . who found that dark matter inside the effective radius is 
negligible with respect to the luminous matter. 

Define p as the ratio of the tidal force and self-gravity. 


IdI IVOtI ^ GM(L)IH aCr-j M(L) 

Ifgravl Ifgravl GnpaRo dr npoH’ 


which reflects the significance of tidal force compared with self¬ 
gravity. In the case of cluster A2744, the mass distribution can 
be modeled by a N EW profile with the parameter r^ ~ 160 kpc 
dZhang et alj|2006l) . We take r 2 Qo, i.e., the radius within which the 
mean density should be 200 xp^, as the boundary of the cluster, and 
set pg = 10 ®Mo/kpc®, the same order of magnitude as the average 
density of the Milky Way. We find that the galaxies with P > I only 
compose 0.1% the volume of the cluster. Thus for a rich cluster like 
A2744, the condition /3 « 1 is satisfied for most galaxies, which 
means that the tidal force is weak compared with the self-gravity 
of th ese galaxies. _ 

ICarter & Luminell dl983ti calculated the response of a star un¬ 
der the tidal force which is weak compared with its self-gravity. 
Analogous to a star, a galaxy will have a trace-free quadrupole 
oscillation, with a frequency Wosc - V^Gpo. Therefore the de¬ 
formation time scale of a galaxy should be Tp =: l/tu^sc = 
1/ VGtrpo. Equation 0 is used to substituting trpo, we obtain 
td == 1 / ^jGM(L)lpH = VStk, where tk^I/^GM(L)IH is the 
Keplerian time scale characterizing the rate of change of the tidal 
force. Since /I <K 1, we have Tq <k Tr. It follows that the galaxy will 
remain close to the stationary equilibriu m state determined by th e 
instantaneous value of the tidal potential dCarter & Lumineil983h . 
where the equilibrium shape of the galaxy will coincide with the 
iso-potential surface of the equilibrium state Oetf = cfigrav + ®x- 
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Figure 1. Upper panel: Illustration of radial orientation angle. Rq = Lq/Ts is 
the dimensionless offset distance, where Lq is the distance between the real 
mass center of the cluster and the reference center used in observation, <p is 
the angle between the direction of the projected elliptical galaxy’s long axis 
and the vector pointing from the galaxy center to the reference center Lower 
panel: Definition of the coordinate system. is the unit vector pointing to 
the direction of the center of cluster, with an angle between which and 
the z axis. I, is the distance from the mass element to vector h is the 
projection of the mass element’s position vector onto Vc 


3 QUANTITATIVE CALCULATION OL THE ELLECT OE 
TIDAL EORCE ON RAD OE GALAXIES IN A 
CLUSTER 

3.1 Coordinate system 

Undoubtedly, an elliptical galaxy with a central luminous ellip¬ 
soid and dark matter halo component is a very complex entity, and 
therefore, we need to assume some necessary simplifications and 
assumptions in order to be able to study mathematically the defor¬ 
mation of such a complicated system under the tidal force from a 
cluster. Since we only consider the tidal distorsion within the effec¬ 
tive radius r ~ 10 kpc, therefore hereafter the effect of dark matter 
is neglected (see section 2 for details). We ass ume that the surface 
of the luminous matter is an oblate ellipsoid l lLambas et al.lll99^ 
at the effective radius, with ellipticity e before tidal distortion, the 


equation of the shape can be written as, 

2 , 2 

+ = 1 


(9) 


where x, y and z axes are the principal axes of the ellipsoid, and v 
and /j are related by the following equation. 


= (1 




( 10 ) 


We require that the vector pointing to the center of the cluster falls 
in the plane x- z, thus the coordinate system can be uniquely de¬ 
fined, as long as the center of the cluster is not on the z axis. 

The direction of the cluster center is denoted as: 


14 = (sin 6c, 0, cos 0c), (H) 

where 9^ is the angle between the direction of the cluster center and 
the z axis. The coordinate system and vector \4 are illustrated in 
the lower panel of Fig.[T] 


3.2 Deformation of a single elliptical galaxy 

In ord er to simplify the calculatio n, according to the homoeoid the¬ 
orem dSinnev & Tremain3ll987h . we assume that the iso-potential 
surface of the gravity potential of the luminous matter at the effec¬ 
tive radius coincides with the local surface of the luminous ellip¬ 
soid. Then the self-gravity potential of an elliptical galaxy at the 
effective radius before tidal distortion can be written as, 

Oi.cav=Al(x"+^I")+A3^^ ( 12 ) 

where Ai = (1 - It is easy to find A 3 =: GnpQ. 

The contribution of tidal force is equivalent to adding a tidal 
potential Ox upon the self-gravity potential Ogca,; here we take, 

= ®grav + ®T- (13) 

A little calculation gives. 

Ox = “■^((1 “ 3cos26c)x^ - 2y^ 

-l-(l -I- 3 cos 26c)z^ + 6xz sin 26c), (14) 

where B = Hence an explicit expression of Ocit is given by: 

®eff = + by^ + cz^ + dxy + exz + fyz, (15) 

where a = Ai - |(1 - 3 cos26c), b = Ai + ^, c = A^ - |(1 + 
3 cos 26c), = “Y sin 26c and d = f = 0. Therefore, the final 
shape of the distorted galaxy is fully described, given the initial e, 
l4, Pg, together wit h the cluster mass pr ofile M(L), which is taken 
as the NFW profile ^Navarro et alJl996h . 

M{r) = 4:?rpor,\ln(l -H r/r^) - ), (16) 

1 -I- r/xs 

where po and the “scale radius” Xj are parameters of the distribu¬ 
tion. Define dg = ^, wherepc is the critical density defined as 
where H is the Hubble constant and G is the gravitational constant, 

^ 200 

do =-, (17) 

3 ln(l-Hc)-c/(l-bc) 

where c is the so-called “concentration parameter” defined as c = 
where X200 is the radius within which the mean density should 
be 200 X Pc ^Navarro et al. | [T^ . Substituting Equation (ED into 
the expression of B, we get 

B = 4;rGpox'^^(ln(l -I- x')----^), (18) 

1 -b X' 

where x' = x/Xj is the dimensionless radius. Given c and x', B can 
be obtained from Equations ( 1171 and ( 1181 . Analogously, Ai, A 3 and 
all the parameters from a to / of the ellipsoid can be calculated. 
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Figure 2. The intrinsic (dashed line) and distorted (solid line) ellipticity 
distributions of galaxies. 


3.3 Projection effects 

From the view of an observer, the ellipsoidal shape of a galaxy is 
projected to become an ellipse. The next step is to determine the 
major axis of the projected ellipse and evaluate the angle between 
that and the direction toward the cluster center. The observer is as¬ 
sumed in the direction Vo, which can be expressed as, 

= (sin 9o cos sin do sin (po, cos 9o). (19) 

The problem is then reduced to such a constrained extremum 
value problem: To find the farthest point i? from the vector ^o in 
the surface of 


Const = aZ + by~ + cZ + dxy + exz +fyz, (20) 

then the projection of i? on the celestial sphere, ^p, is the direction 
of the major axis. Define V^p as the projection of vector on the 
celestial sphere, we obtain the radial orientation angle as. 


if = arccos 


I i^cp II4 r 


( 21 ) 


The details of solving the constrained extremum value problem are 
presented in the Appendix. 


3.4 RAD of elliptical galaxies 

In order to derive RAD of elliptical galaxies in a cluster, we use 
Monte-Carlo simulations. In each run, a set of elliptical galaxies 
are generated, each of which is given a set of parameters sam¬ 
pled with certain probability distributions, as discussed below. With 
M{r) given, each galaxy from the set results in a i/j, {i = 1,2...A), 
c = 10 is used to simulate a cluster with high concentration like 
Abell 2744 (see secti on 4.1); note that c = 10 is also the concentra¬ 
tion of Coma cluster l lRines & Diaferiol2006l) where RA was found 
for the first time. The effective density of each galaxy is set to be 
10* times po- In this run, 100,000 galaxies are generated, with an 
intrinsic ellipticity distrib ution in Fig.[3 which is the ellipticity dis¬ 
tribution of field galaxies l lFasano & Violl99ll) . Here the ellipticity 
of a generated ellipsoid is defined as e = , where /max and 

/min are the major and minor axes of the ellipsoid, respectively. The 
dimensionless projected distance of galaxies are sampled with a 
probability inversely proportional to r’^, from = 2.5 to = 5, 




Radial Angle (”) 


Figure 3. RADs for different center offset distances Rq. 


corresponding to an appro ximation to th e projected galaxies num¬ 
ber density distribution in iBaierl (Il976l) . Under such parameters, 
the maximum value of fj is 0.0057, which is far smaller than unity. 
For each galaxy, the direction toward the cluster center Vc and the 
direction toward the observer Vo in the principal axes coordinate 
system are sampled with an isotropic probability, corresponding to 
an isotropic position distribution, which is assumed as the initial 
configuration of galaxies in the cluster. 

Fig.[3]shows RADs for the different center offset distances Rq, 
where Ro = To/A is defined as the dimensionless offset distance, 
where Lq is the distance between the real mass center of the clus¬ 
ter and the reference center used in observation. When Rq = 0, 
the upper-left panel in Fig.[3shows a single peak. Therefore if the 
reference center is the real cluster center, the major axes of the el¬ 
liptical galaxies tend to point to the cluster center, i.e., we can see 
obvious radial alignment of the elliptical galaxies. However if the 
reference center has an offset distance from the real cluster center, 
since the major axes tend to point to the real center rather than the 
reference center, two peaks are found in RADs. As the center devia¬ 
tion distance increases, the two peaks depart from each other farther 
with increasing Rq. Then the two peaks merge at -90/90 degrees 
as Ro keeps increasing. Finally, RAD becomes uniform distribution 
if Rq s> 0 , and at this moment, radial angle distribution becomes 
position angle distribution. 

We also plot the distorted ellipticity distribution of ellipsoids 
after the simulation in Fig.|2] The width of the ellipticity distribu¬ 
tion is broadened about Aw = W 2 -wi 0.025, where wi and W 2 are 
the full widths at the half maximum (FWHM) of the ellipticity dis¬ 
tributions before and after the simulation, respectively; meanwhile, 
the ellipticity distribution is shifted to higher values compared with 
the intrinsic one. 

In order to further reveal the effect when the center of the clus¬ 
ter is mis-identified observationally, we define the “amplitude” of 
RAD as Arad = , where Amax, Amin and A are the highest 

point, lowest point and mean value of RAD. Arad intuitively sug¬ 
gests the significance that a peak-type distribution deviates from 
uniform distribution, i.e., Arad - 0 when the distribution is simi- 
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Ro 


Figure 4. Amplitude of RAD versus Rq. 

lar to uniform distribution, while » 0 if there are significant 
peaks in the distribution. Fig.|4]shows A^ad as a function of the offset 
distance from the reference center to the true center of the cluster 
mass distribution. On the whole, A^ad decreases with increasing Rg. 
For Rq » 1, a radial angle is equivalent to a position angle. The 
result in Fig. |4] explains why it is difficult to find RA with only the 
position angle distrihution of elliptical galaxies in a cluster. 

4 OBSERVED RAD OF ABELL 2744 
4.1 Selection of Elliptical Galaxies 

Ahell 2744 is one of the richest clusters of galaxies at intermedi¬ 
ate re dshift z ~ 0.308, known as a gravitational lens (lMerten_e£^ 
20^. Its parameters are estimated as - 16 0 knc l lZhang et al. 
20061) . r 2 oo = 2.506 Mpc jPemarco et alj|2003l) . hence the concen¬ 
tration parameter c = r 2 oo/A - 15.7. Such a high mass concentra¬ 
tion is expected to show some signature in its RAD. We thus choose 
this cluster to make the comparison with our theoretical work. 
The coordinate system is cent ered at a (J2000.0) = OO''14™21.04% 
6 (J2000.0) = -30°23“52.4’^ l lBoschin et al.ll2012t). However, not e 
that this cluster is an actively merging cluster dOwers et alJl201 ih : 
therefore, the cluster central position cannot he determined unam¬ 
biguously l lZhang et al.ll2006l) . which may bias RAD. 

We retrieve HST (ACS/WFC) images in the F606W and 
F814W hands of two fields in the immediate vicinities of the clus¬ 
ter center (Date: 10/27/2009, PID; 11689, PI: Dupke). The asso¬ 
ciated images of A2744 are fully processed and drizzled. Pho- 
tometry for ohjects is c arried out using the Sextractor package 
dBertin & Amoutsll 19961) . The configuration parameters are listed 
in Tahle[T] and the magnitudes of the sources in the cluster are cal¬ 
ibrated by the MAG_ZEROPOINT parameter in the AB system. 
The derived source catalogs in the two bands are matched to obtain 
3005 true objects. The parameter Class/Star> 0.9 is used to get rid 
of stars, and then the removed stars are visually inspected in the im¬ 
ages. The arcs and arclets, the spurious objects, the sources in the 
margins of images, and the sources inside other giant bright sources 
are visually identified and removed. The optical magnitudes of tar¬ 
gets are corrected for the foreground extinction from the Galaxy ac- 
cording to the Schlegel-Fink beiner-Davis Galactic reddening map 
dSchlaflv & FinkbeinejbOlII) . 


Parameters Values 


DETECT J4INAREA 

32 

DETECT.THRESH 

1.3 

DEBLEND JSTHRESH 

64 

DEBLEND JAINCONT 

0.005 

CLEANJARAM 

1.2 

BACK.SIZE 

40 

BACK.FILTERSIZE 

3 

BACKPHOTO.THICK 

24 


Table 1. Sextractor configuration parameters for Abell 2744. 
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Figure 5. CMD of Abell 2744. The galaxies with color -0.037m8i4 4- 

1.610/q 2 m and m 8 i 4 < 23.0 m are considered as the E/SO galaxies and 
denoted by a box (red in the online aificle). 


Limited by the ACS field size, the measurement reaches out 
up to a distance of 1.28 Mpc from the cluster center, which is about 
eight times of r^. These objects include both the foreground and 
background galaxies. In order to select the real cluster members, a 
color-magnitude diagram (CMD) is plotted in Fig. [fusing the mag¬ 
nitudes in the two filter bands, i.e., color = mag = mgoe- 

Since the colors of elliptical galaxies are redder than spiral galaxies 
jBaldrv & Glazebrookll2005[) . and the Butcher-Oemler condition is 
satisfied, i.e., the blu e galaxies are those at lea st 0.2 mag bluer than 
the cluster ridgeline jButcher & Oemle^ll984h . a linear fitting, i.e., 

color - k ■ mag 4- d, (22) 

where k, d are the slope and interce pt respectively, is app lied to 
obtain the slope of the red sequence l lBrow et al.lll992arbh in the 
CMD, which is known as the collection of E/SO galaxies. Here 
k ^ -0.037, d 1.610 are obtained; therefore we select the E/SO 

galaxies in the color range of color -0.037/n8i4 4- l.hlO/gj m. A 

faint end at about ra 8 i 4 ~ 23.0 m is used to di vide the member ellip - 
tical galaxies from the background galaxies dRomano et al ]|2 oI 3), 
which is approximately consistent with the rest-frame r-band ab¬ 
solute magnitude of ~ -17 m. Einally 199 elliptical galaxies 
are obtained. These selected elliptical galaxies are then visually in¬ 
spected in the images. 
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Figure 6. Observed RADs of Abell 2744. The upper panel is RAD of the 
all sources; the bottom one is RAD of the elliptical galaxies. RAD of the el¬ 
liptical galaxies is fitted by a double-Gaussian component plus a Lorentzian 
component. 


4.2 RAD of Abell 2744 

With SExtractor, we are able to quantitatively measure the orienta¬ 
tions of the detected objects. We define the radial angle ^ = 0 when 
the major axis of a projected ellipsoid points to the cluster center, 
and (fi = ±90° when the major axis orientates along the tangential 
direction. By separating the range of -90° —1-90° into 18 bins, we 
count the number of members in each bin. In order to obtain RAD 
with a high significance, we add together the number of galaxies 
in each bin of ip from the two filter bands. RAD of the all sources 
of Abell 2744, including the foreground and background galaxies, 
and RAD of the elliptical galaxies in Abell 2744 are presented in 
Fig.|6] 

There are apparent tangential peaks centering at about ±90° 
in the upper panel of Fig. which should be produced by grav¬ 
itational lensing (e.g. ISinailet ^1 19971 ; iMerten et al.ll201lh . since 
the shap es of backgrou nd galaxies are stretched tangentially after 
lensing djoachimll 19971) . Further more, since the selected elliptical 
galaxies are contaminated by some background galaxies, RAD in 
the bottom panel of Fig.|^also show weak tangential peaks, which 
are fitted by a double-Gaussian component, 

VZtrcTG Vt^ctg 

where pQ = -90° is fixed. Additionally, a peak is located at close 


iVL 

8.3(1.7)x 10- 

0.17(4.11) 

CM 

DG-fL 

Ng 

3.0(1.0)x 10- 

‘PG 

-90(fixed) 

o-G 

22(8) 

X^ldof 

14.0/13 




DG+C 





Nq 


Ng 

PG 

o-G 

yidof 


4.5(0.6) 


5.3(5.4) X 10 

-90(fixed) 

6.8(9.5) 

27.8/15 


Table 2. Fitting results of RAD of the elliptical galaxies. DG denotes a 
double-Gaussian component, L denotes a Lorentzian component, and C de¬ 
notes a constant component. All errors are at Icr throughout this work. 


to ^ = 0, which is fitted by a Lorentzian component, 

Ai <Tr 

^=^7 - 2T (24) 

n[(ip - ifi^y + o-L^] 

As a comparison, we also fit the central peak with a constant, 

A = Ac. (25) 

The two different fitting results are listed in Table |2] The rela¬ 
tive error of Ac is 20%, suggesting that the central peak really 
exists, and thus the orientations of the elliptical galaxies are not 
random. F-test is performed to obtain the probability p that the 
model DG±L is better than the model DG±C. Given the value of 
F,F =xr doftlxl ■ ■where dofi = 13, rfo /2 = 15, mdx\,xl 
are the variances of the two models, respectively, thus we obtain 
p ~ 83%. 


5 CONCLUSION AND DISCUSSION 

We estimated qualitatively the shape changing time scale of a 
galaxy in a cluster, which is found to be much shorter than its Kep- 
lerian time scale around the center of the cluster for most member 
galaxies. Thus we ignored the dynamical process of shape chang¬ 
ing and assumed an equilibrium shape of galaxies. The shape of 
an elliptical galaxy under the tidal force of a cluster is calculated, 
by assuming that the density contour of the galaxy coincides with 
the equipotential surface of the sum of the tidal potential and self- 
gravitational potential of the galaxy. We then used the Monte-Carlo 
method to simulate RAD of galaxies in a cluster, taking the NFW 
mass profile of the cluster and the initial ellipticity distribution as 
that of field galaxies. The Monte-Carlo simulations also find that 
the elliptical galaxies in a cluster tend to be radially aligned if the 
real physical center is chosen as the reference center. If the ref¬ 
erence center has an offset distance from the real center, then the 
farther the reference point away from the cluster center, the weaker 
the anisotropy of RAD will be shown. Therefore in order to find 
significant RA, the reference center should be set to the real phys¬ 
ical center of the cluster. In addition, the distribution of position 
angles which is equivalent to Ro » 0 cannot reflect RA. The ob¬ 
served RAD of cluster Abell 2744 is presented, using the data from 
HST. Comparing the observed RADs in Fig. with the simulated 
ones in Fig. [3 we conclude that the observed non-uniform RADs 
are due to the tidal effects of these clusters. 

In our model, the gravitational interaction between two nearby 
galaxies is neglected. We treated the galaxies as if they were mov¬ 
ing under a smooth averaged field. In fact when two galaxies get 
close, the tidal effects from the nearby galaxy will weight over the 
entire cluster. This will increase the randomness of RAD of the 
whole cluster, because encounter events are random and the ef¬ 
fects are local. The non-spherical symmetry and substructures of 
the cluster are also not included, since we applied the isotropic 
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NFW mass profile. Those effects will be included in our follow¬ 
ing works. Our simulations place galaxies in the region where the 
self-gravity of the galaxies are far greater than the tidal force from 
the cluster, so tidal disruption of galaxies has not been taken into 
consideration. 

In section 2, we assumed that a galaxy has a trace-free 
quadrupole oscillation u nder the weak tidal distorsion like a star 
dCarter & Lumine3ll983h with a period Tosc ~ 0.1 Gyr. However, 
an elliptical galaxy is a stellar system rather than a fluid system, 
and the oscillation can be quickly damped by Landau damping 
dSinnev & Tremain^ll987l : Iweinberdll994li . Therefore, the time 
needed from the beginning of deformation to finally settling on 
the equilibrium state is indeed the damping time scale Tjampi if 
Tdamp > fo.sc- In a finite stellar system like an elli ptical galaxy, the 
time scale of damping the quadrupole oscillation dWeinberdl 19941 
I = 2 mode) is shorter than Tosc/ 0.608 - l.brosc (see Weinber j 
d 19941) for details). Therefore an elliptical galaxy will remain on 
the stationary equilibrium after 1.6 tosc, which is also shorter than 
Tk. Therefore Oefp = ®grav -t Ox is a good approximation, where Op 
is instantaneous. 

In the Monte-Carlo simulations, the initial ellipticity distribu¬ 
tion and the effective density of the galaxies are chosen as those 
of field galaxies. There is evidence showing that differences might 
exist between th e density of cluster galaxies and field galaxies 
dPu & Hanll20TIh : it is however possible that these differences are 
caused by the tidal effect we calculated here, which may change 
the density and ellipticity distributions of the member galaxies of 
a cluster. Although different initial ellipticity distribution and ef¬ 
fective density chosen do not change our conclusion qualitatively, 
these will become important when analyzing RAD quantitatively. 
The effects of weak lensing are also not taken into consideration, 
which would cause images of some of galaxies in the far and back 
side of the cluster stretched along the tangential direction, and thus 
may weaken the observed RA by the tidal force of the cluster. 

Finally, in this work we aim to study the possibility of radial 
alignment of elliptical galaxies, the original orientations of which 
are random, under the tidal force in a cluster. There are certainly 
other possibilities of radial alignment; for example, the galaxies en¬ 
tering the cluster could be radially pre-aligned, since tidal forces are 
ubiquitous. However, the issue is beyond the scope of the present 
work, and could be an interesting topic of future investigations. It 
is possible to even infer the original orientation and ellipticity dis¬ 
tributions of galaxies by studying the observed distributions using 
our analysis method, if the gravitational field of the cluster is de¬ 
termined independently with other methods. Once again, this could 
be an interesting topic of future investigations. 
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Substitute j of Equation l lAlOt into Equation l lA4k we can get both 
the major and minor axes of the projected ellipse. 


APPENDIX A: DETERMINING THE MAJOR AXIS OE AN 
ELLIPSOID PROJECTED ONTO THE CELESTIAL 
PLANE 


The projection of an ellipsoid is an ellipse. The problem is how to 
determine the major axis of the ellipse, given the equation of the 
ellipsoid and the direction along which the ellipsoid is projected. 
An ellipsoid described by equation, 

ax^ + by^ + cZ -I- dxy ■¥ exz + fyz - 1, (Al) 


can be reexpressed in a matrix form, 


where ^ = (x, y, z) and A = 

Suppose the direction of the observer is 


ii-1 

= 0, 


(A2) 

a 

d/2 

e/2 ' 


d/2 

b 

f/2 


e/2 

fl2 

c , 



n = (sindcos sindsini^, cos 6), (A3) 


so the problem is to find a point on the ellipsoid farthest from the 
vector d. 

Define = Ad + pv, which is a bunch of rays that are parallel 
to n, where p and A are real numbers and if is any unit vector that is 
perpendicular to it. v can be written as, 

, (- cos 0 -I- sin 0 sin (b /, - sin 6 cos dj /, sin 9 cos 0) 

i7= ■' - - . (A4) 

Jp sin^ 6 - 2y sin 6 cos d sin 0 -l- cos^ 9 + cos^ 0 sin^ 9 
where j e (-oo, -t-oo). 

The intersection points of the ellipsoid and the rays bunch sat¬ 
isfy the ellipsoid equation, 

-1 = 0. (A5) 

Taken the expression of 14,, Equation l lA5l l can be written as, 

PrPAd + p^iPAiS + lAprPAij- 1=0, (A6) 

which is a quadratic function of A. The ray Ci, is tangential to the 
ellipsoid when Equation l lA6t has multiple roots, i.e., 

= (fAd)-'. (A7) 


The next step is to find j that maximizes p. Using the expression of 
A and v. Equation ( IA7t is rewritten as. 


where 


f + aj + b 
f+aj + F 


a = -2 cot 9 sin 0 
b = cos^ 0 + cot^ 9 



(A8) 


(A9) 


and (7a, Ut, and D above can also be expressed explicitly in term of 
a, b, c, d, e, /, 9, cfi, in a lengthy but trivial way. 











